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ABSTRACT 



Aims. Previous studies of the X-ray flux and spectral properties of IRXS J 170849-4009 10 showed a possible correlation with the 
spin glitches that occurred in 1999 and 2001. However, due to the sparseness of spectral measurements and the paucity of detected 
glitches, no firm conclusion could be drawn. 

Methods. We retrieved and analysed archival Rossi-XTE pointings of IRXS J 170849-4009 10 covering the time interval between 
January 2003 and June 2006 and carried out a detailed timing analysis with phase fitting techniques. 

Results. We detected two large glitches of 1.2 and 2.1x10"^) that occurred in January and June 2005. The occurrence times of 
these glitches are in agreement with the predictions made in our previous studies. This finding strongly suggests a connection between 
the flux, spectral and timing properties of IRXS J170849-400910. 
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1. Introduction 

IRXS J170849-400910 (hereafter RXSJ1708) is one of the 
anomalous X-ray pulsars (AXPs), a small group of peculiar neu- 
tron stars (NSs) that are currently believed to have super-strong 
magnetic fields, B ~ 10'"* - 10'^ G, hence dubbed "magne- 
tars" (though other possible models are no t completely ruled out 
by observations; for a recent review see Woods & Thompsod 
, 2006 ). RXSJ1708 was first discovered by ROSAT ( Voges et aU 
ri996l) . while ~ 1 1 s coheren t pulsations were detected with 
ASCA dSugizaki et al.L Il997h . Early measurements suggested 
that it was a fairly stable rotator, with a spin-down r ate of 

~ 2 X 10^^'ss and a soft spectrum dlsrael et all 1 19991 

Events of sudden spin-up (glitches) with very different 
post-glitch recovery were detect e d in RXSJ1708 by Rossi-XTE 
in 1999 and 2001 (|l£asBieLaU l2000t iDalFOsso et all l200l 
iKaspi & Gavrii]ll2003h . The rather short interglitch time makes 
this AXP a fr equent glitcher am ong neutron stars. 

Recently. iRea et afl (l2005ah noticed that the long term (over 
5 years) variations in the source X-ray flux and spectral hardness 
are correlated, with both quantities reaching a maximum close 
to the epochs of the two glitches in 1999 and 2001. The addition 
of new data obtained from Chandra (2004) an d Swift observa- 
tions confirmed the flux-hardening correlation ( Campana et aLl 
l2007h . extending its val idity to har d X-ra ys jGotz et all l2007l) . 
Following these studies. IZane et al.l (|2007|) and Rea et a l1 (l2007h 
suggested that the long-term variations may have a cyclic be- 
havior with a recurrence time of x 5-lOyrs, possibly due 
to a periodic twisting/untwisting of the star magnetosphere 
jBeloborodov & Thompson! |2007[) . Correspondingly, the source 
was expected to re-enter into a glitching active phase during 
2005-2006, close to the latest maximum in the source flux. 
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Thus, we analyzed the Rossi-XTE archival data spanning 
the latest 3.5 years and performed a phase coherent analysis 
of the pulse arrival times in a search for new glitches. We de- 
tected two large glitches occurring on J anuary 2005 and May 
2005, in agreem ent with the prediction of IZane et al.l (l2007h and 
iRea et al] (I2007D and correlating with the flux and spectral evo- 
lution. 



2. Rossi-XTE observations and timing analysis 

We analyzed 204 Rossi-XTE archival observations of 
RXSJ 170^1 These span from 2003 January 5th to 2006 
June 3rd. We restri cted our analysis to the PCA instrument 
(iJahoda et al.L 11996 ^ which was operated in good Xenon data 
mode with a time resolution of 1 fis and 256 energy bins between 
2 and 120 keV. Raw data were reduced using the ftools v6.2, 
provided by the High Energy Astrophysics Science Archive 
Research Centefl The events were extracted in the 2.5-16 keV 
energy range, and binned into light curves of 0.125 s resolution. 
Photon arrival time correction to the barycentric dynamical time 
(TDB) was applied by usi ng fxbary and th e (~1" accurate) 
source position provided by llsrael et al.l (l2003h . 

A phase-coherent timing solution was first obtained using 
a relatively long (29ks) archival Chandra observations caiTied 
out on 2004 July 3rd (starting time at MJD 53189.151 489 TDB; 
for d etails on the data reduction of this dataset see iRea et al.L 
l2005h) . This provided a sufliciently accurate period determi- 
nation, P= 11.00231(3)s, such that no pulse cycle was missed 
when extrapolating this value to the epoch of the closest Rossi- 
XTE pointing (2004 July 1st). A phase-coherent timing solu- 
tion was inferred in the time interval between 2004 May 1st 

' ftp://legacy.gsfc.nasa.gov/xte/data/archive 

^ |http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftooIs_menu.html| 
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Table 1: Measured parameters for the two glitches detected in 
the Rossi-XTE data of RXSJ1708. Icr errors in the last digit are 
quoted in parenthesis. 
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Fig. 1: Rossi-XTE time residuals for the time interval January 
2003 - June 2006 after subtraction of the phase-coherent P-P 
solution inferred bv lDair Osso et al ] 12003). The inset shows the 
time interval over which we detected the two glitches. 
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Fig. 2: Time residuals of Rossi-XTE RXSJ1708 observations 
from June 2004 to June 2006 after subtraction of pre- and 
post-glitch model N.l and the polynomial post-glitch model 
N.2. Vertical hnes indicate the assumed/inferred epochs of the 
gUtches. 



and November 16th corresponding to v= 0.090890035(1) Hz and 
v= -1.5884(14) (epoch 53819.0 MJD; Icrc.l. are reported). 

The inclusion of the 2005 and 2006 datasets showed large 
disagreement with the above inferred timing solution, with two 
evident "jumps" in phase both marking a decrease in the pe- 
riod value (see blue and dark violet filled circles in Figure 
[T]l, strongly suggesting that two further glitches occurred at 
the end of 2004 / begin 2005 (MJD 53370) and in May 2005 
(MJD 53550). No signature for similarly large glitches was in- 
stead found in the datasets spanning between January 2003 and 
Jun e 2004. The phases rep orted in Figure [T] are obtained using 
the iDairOsso et al.l (120031) 2001 post-gHtch solution (horizon- 
tal solid line) which is nearly coincident to the phase coher- 
ent solution inferred above. Thus, we applied a detailed tim- 
ing analysis to phase residuals were the two large glitches oc- 
curred. A more exhaustive study and modelling of the whole 
phase and glitch history of RXSJ1708 will be presented else- 
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where ("Dair Osso et al.', '2007"). Following the analysis scheme 
outlined in Dall'Osso et al. (2003) we inferred the main param- 
eters of the two detected large glitches (see Table [U. Figure |2] 
shows the phase residuals after subtraction of the timing solu- 
tions pre- and post-glitch N.l, and post-glitch N.2. Both the 
newly detected glitches reveal large jumps in the spindown rate, 
Av/v ~ 1 X 10"^ and Av/v ~ -0.1, among the largest ever 
observed in glitches that lack a significant short-term exponen- 
tial recovery. Remarkably, they have opposite signs: the second 
glitch has cancelled the effect of the previous increase in v and, 
actually, somewhat overshot it. 

The jump in spin frequency after the first glitch appears have 
been recovered in ~ 120 d. The upper limit on v after the first 
glitch implies that the jump in v could have been recovered, if 
at all, only on a much longer timescale and not until the second 
glitch occurred, ~ 175 d after the first one. At the second glitch, 
an even larger spin-up occurred, accompanied by a significant 
flattening of the spindown trend. Thus, the spin up started with a 
sudden increase and then it slowly continued. A large and nega- 
tive second derivative is found to be highly significant in the fits, 
which corresponds to a long-term recovery of Av on a timescale 
^ (503 + 48) d. We note that the v term detected after the second 
glitch brings v back to the value it had before the first glitch in 
~ (159 + 26) d, and to steeper values afterwards. The extrapo- 
lated spin frequency at 500 d from the second glitch gives an ad- 
ditional spin-up ~ (5.7 + 0.5) x 10"^ s"', while at the same epoch 
V returned to the value v(Af = 500d) ^ -(1.70 ± 0.02) x 10"'^ 
s"^ that it had after the first glitch. 

3. Discussion and conclusions 

RXSJ1708 experienced two new glitches, both with a large frac- 
tional amplitude, Av/v ~ 1.2 x 10"* and 2.1 x 10"* respectively, 
comparable to the so-called giant glitches observed from Vela 
and to the gl itch previ ously detected from thi s sourc e in May 
2001 (Dafl'O sso et"an. 2003; Kaspi & Gavriill l2003h . This re- 
sult strongly suggests that giant glitches seem to be the rule for 
this source. 

It has been noticed (see iDalFOsso et al.L 120031; IWang et all 

2000, for details) that most glitch models are diflicult to recon- 
cile with observations of a growing number of glitches. In par- 
ticular, all models requiring a catastrophic (i.e. widespread) un- 
pin ning of crus t al superfluid (such as the vortex creep models, 
see lAbar et al.L IT9841 11993L Il989h seem less promising, since 
they require simple correlations between the amplitudes and re- 
covery timescales of glitches in a single source as well as an 
approxi mately constant Q-value (where Q is the recovery frac- 
tion, see iDairOsso et al.l (l2003l) ). which is found not to be the 
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case, in general, for radiopulsars. Although these problems can 
be partially solved by invoking local unpinning, this appears to 
be a rather ad hoc assumption (Jones, 2002). Furthermore, it has 
been recently pointed out that the observation of long-period pre- 
cession i n a few puls ars is incompatible with pinning of crustal 
vortices (lLinkll2006l) . 

Models where the internal angular momentum reservoir is 
in the NS core have some advantages with respect to the above 
points and can in principle explain a wider range of glitch prop- 
erties (e.g.see lDairOsso et alll2003i for a more complete dis- 
cussion). In this respect, it is worth noticing the unusually large 
jumps in the spindown rate of the two new events. A ~ 7% in- 
crease found after the first glitch would rule out the crustal su- 
perfluid as a momentum reservoir: the core superfluid - and a 
remarkable fraction of it - would necessarily be required to have 
been involved in the glitch. On the other hand, the large flat- 
tening that follows the second glitch is also very peculiar (al- 
though a residual steepening of the spindown trend, not recov- 
ered promptly, wa s also found in the May 2001 glitch from this 
same source, Da lFOsso et all l2003h . In general, it is not easy 
to explain long-term ofi'sets from the secular spindown trend 
in ghtch m odels based solely on crustal superfluid instabilities 
(lLinketal.L[T992l) . since these require permanent decoupling of 
some fraction of the internal superfluid from the bulk of the NS 
mass. More extreme difliculties are posed by the slowly recov- 
ered decrease in the spindown rate found in the second glitch: 
this apparently implies either a re-coupling of an internal com- 
ponent that was previously decoupled from the crust, or a corre- 
sponding decrease in the spindown torque. 

A potential argument against core-based models arises when 
using Vgi (the average spin-up rate caused by glitches) to es- 
timate the fractional mome nt of inerti a of the in ternal reser- 
voir o f angular momentum dLvne et akl [2000; Dal FOsso et al.L 
l2003h . The parameters of the three large glitches detected so far 
from RXSJ1708 yield Vj, ^ 2.65 x 10"'^ s^^. Since the ratio 
Vgilv < Ires/Ic, onc dcrivcs Ireslh ^ 0.017 whcrc Ires and Ic 
are the reservoir and stellar moments of inertia. This is compa- 
rable with the value inferred typically in glitching radio pulsars, 
and consi stent with the cru stal superfluid being the momentum 
reservoir dLvne et aill2000l) . However, we notice that the equal- 
ity in the above formula obtains in the ideal case where the in- 
ternal reservoir is completely decoupled from the crust. In the 
presence of coupling this component would also spin-down in 
between glitches at a rate Vyes < v. Therefore, the following 
equality, Vgi/{v - v^i) = heslh should be considered instead. If 
glitches are produced once a critical rotational lag Vcr is reached, 
then V - Vres ~ Vcr/tg, whcrc tg is the typical interglitch time. 
If we assume I,-es - 0.1/c, as required to explain the changes in 
V found for the two new glitches, the measured Vgi and a typical 
interglitch time <2ys, we find Vcr < 1.6x 10"^ s ' . Whether such 
a low critical lag is meani ngful on physical grou nds is matter for 
future investigation (see lDalFOsso et £11120031 for similar con- 
clusions). An alternative explanation, in terms of starquakes and 
subsequent movements of cracked platelets, has been proposed 
for at least the large glitches of radio pulsars ( iRudermanl 119911: 
[Ruderman et al., 1998). Quakes may be caused by the stress act- 
ing on the crust following the interaction in the core between 
superfluid neutrons and magnetic flux tubes threading the crust. 
Alternatively, the crustal strain that is accumulated during the 
growing of a twist in a magnetar's magnetosphere may repre- 
sent a trigger in AXPs. Only the moving sector and its surround- 
ings are affec ted, leaving the rest of the star mainly unperturbed 
dJonesl I2OO2I) . The "local nature" of starquake models make 



> 

CD 



4.5 
4 

2 3.5 
I 

^ 3 
3.2 



X 
CD 

T3 2.8 



o 2.6 
o 



2.4 - 



2.2 



H h 



I I 



H 1 1 h 



-l-H- 



H — h 



5.2) 



:x104 5.3x104 5.4x104 
Time (MJD) 

Fig . 3 : Long-term spectral history of RXS J 1708 showing the cor- 
related hardness-intensity variation. All reported fluxes are ab- 
sorbed (with Nh - 1.36 X 10^ ^ cm~^ for a ll the spectra) and in 
the 1-lOkeV energy range (see lGotz et al, l(2007) for further de- 
tails). Vertical dashed lines mark the glitches detected in 1999, 
2001 , January and May 2005, from left to right respectively. Flux 
cross calibrations among instruments are less than 15% (note 
that the latest two measurements are both inferred from Swift 
XRT and implies a flux variability of more than 25%. 



them promising for re producing the large var iety of properties 
found in AXP g litches d PaU'Osso et an i2003l). 

In lRea et"al] (l2005al) . we proposed that the observed correla- 
tion between the X-ray flux and the spectral hardness may be ex- 
plained if the evolution is regulated by the cha nge of a "twist" in 
the m agnetosphere (see also Th ompson. Lvutikov & Kulkarnil 
i2002h . The evolving magnetic field is expected to fracture the 
crust at intervals, eventually causing an increased activity and 
large amplitude glitches. We found that observations collected 
until 2003 were consistent with a scenario in which the twist an- 
gle was steadily increasing before the glitch epochs, culminating 
with glitches and a period of increased timing noise, and then 
decreasing, leading to a smaller flux and a softer spectrum (see 
Figure[3]l. Interestingly, the same model provides a natural expla- 
nation for the new period of glitching activity, that was foreseen 
in our previous papers. Nonetheless, we emphasise that while we 
do expect glitching activity corresponding of an increasing stress 
of the crust caused by a growing twist, glitches might also occur 
outside these epochs, in particular if, as stressed before, glitches 
with different properties (such as amplitude and recovery) may 
reflect a difference in triggering mechanism. We tried to esti- 
mate a reliable false alarm probability for the observed glitches 
to occur by chance close to the flux maxima displayed by ob- 
servations, by assuming a uniform probability of flux values at 
glitch epochs. The false alarm probability depends on the max- 
imum assumed flux range (2.5xl0""erg cm"^ s"', 60% more 
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than that observed; see Figure [3j and the given flux threshold Zane, S. et al. 2007, in "The multicoloured landscape of compact objects and 
value (4.2xlO"''ergcm"^S"') above which glitches OCCufl The their explosive origin", Cefalu' (Sicily), 11-24 June 2006, AIP, in press 

above, somewhat conservative assumptions give a false alarm 
probability of 3%. Based on the above considerations we believe 
that the proposed link between timing and spectral properties of 
RXSJ1708 is rather good. However, given the poor knowledge 
of the true flux history and the not statistically independent na- 
ture of glitches, we believe that only future observations will be 
able to unambiguously confirm the link suggested by observa- 
tions. Nonetheless, when the 1999 glitch is also added to the 
above considerations, the false alarm probability decreases to 
1%. Thus, if such correlated long-term variation will be further 
confirmed and discovered in other AXPs/SGRs, its X-ray moni- 
toring might become an excellent tool to gather a deeper under- 
standing of magnetars and , more gen erally, neutron star glitches. 

Finally, we note that iDib et al. I (2007) submitted a paper 
on an independent discovery of glitches from RXSJ1708 in the 
same dataset (after our original report at the Seattle Workshop on 
27th June 200'fl) : they considered our glitch N. 1 as a candidate 
rather than a true glitch. A detailed comparison between their 
results and ou rs is indeed premature and beyond the aims of the 
present work dPalFOsso et al.Ll2007l) . None theless, we note t hat 
a source of discrepancy is due to the use bv lDib et al. I (l2007l) of 
high-order frequency derivatives (used to identify glitches) in the 
presence of gaps in the phase-series. It is indeed not surprising 
that tw o out of the three candidate glitches reported bv .Dib et alTI 
(l2007h are found when data gaps are also present. 
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